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Studies of 5-singuIar vortex dynamics on a rotating sphere are relatively rare (see [1,2] and 
references in [2]) though they can be useful for modeling of different hydrodynamic, geophysical and 
astrophysical systems, in particular for modeling of atmospheric centers of action [3, 4] and blockings 
[5-7]. It should be noted here a problem of infinite dimension of usually analyzed dynamic systems. 
Such a problem appears as a result of the invariant absolute vorticity representation (in [1, 2, 6] 
particularly) on a rotating sphere [8] as a formal superposition of 5-singular point vortices (local 
vorticity) and background global vorticity (associated with the Earth rotation). Interaction between 
local and global vorticity (due to the hydrodynamics equations nonlinearity) leads to the structural 
noninvariance of the local vortices and to the necessity of the approximate considerations. 

In this paper, strong interaction of the local and global vorticity due to the introduction of 5- 
singular (localized) representation of the invariant absolute vorticity is exactly considered. It 
eliminates the cause of the infinite dimension problem for the appropriate dynamic system. This 
representation allows to get a new exact weak solution for an ideal incompressible fiuid 
hydrodynamics on a rotating sphere. Such a solution has a form of 2N-dimensional Hamiltonian 
dynamic system for arbitrary number N>2 of strongly interacting point diametrically conjugated 
vortex pairs (CVP) consisting of the point vortexes with the same absolute value of vorticity and 
different circulation signs that are placed in diametrically conjugated sphere points. It is shown that 
one such CVP is the only one allowed by hydrodynamic equations elementary S -singular stationary 
vortex object on a sphere unlike the sole 5 -singular vortex object introduced in [1]. 

Existence for S -singular vortices of a stable absolute stationary pair (N=2) vortical 
cyclonic-anticyclonic solution (not considered earlier) was determined. Such solution for the vortical 
system in the Northern Hemisphere is possible only in the case with the more intense and more high- 
latitude cyclonic vortex than the anticyclonic vortex at the same longitude. It agrees with observations 
for coupled atmospheric centers of action (ACA) in the Northern Hemisphere [3]. Possibility of 

modeling of the splitting-type blocking [5] with a more intensive anticyclone to the north with 
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respect to the cyclone is also considered. Contrary to the absolute stationary mode of the coupled 
vortex ACA systems the only relative stationary mode is realized with a finite velocity motion to the 
west of vortex pair characterizing blocking (with significantly less intensive vortices than for ACA). 
Significant influence of polar vortices on the stability of the coupled vortex ACA system and vortex 
pairs simulating blocking is found. Comparison of theoretical results with observed and reanalysis 
data [3] for coupled quasi-stationary ACA systems over the oceans in the NH is conducted. Exact 
solution obtained for N=2 is used for analysis of approximate results [6] in relation to the vortex 
modeling problem of atmospheric blockings. 



1. From three-dimensional hydrodynamic equations for rotating with angular velocity Q=const ideal 
incompressible fluid in the absence of radial fluxes with Fr=0 it follows absolute vorticity invariance 
(0 = (0^ -i-2Qcos6' [8]: 

Deo _dco ^0 dco ^(p dco _r. ..^ 
Dt 8t r 80 rsmOdcp 

Here (r, 6, a) are the spherical coordinates, Vfo =- — = — . ^ are components of the non- 

^ r ^ rsmO dq) 

divergent velocity field, y/ is the stream function related to cor, and cor is the radial local vorticity 
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component =-A^^/^,A s — ( — sisP — + y) is the spherical part of the Laplace 



operator. 

Let find exact weak solution of (1) for the following 6-singular representation of the absolute 
vorticity field co on a sphere of radius r=R in the form: 

i?^ i=lR'^sm0i 

where r^- are intensities of CVP, characterized by the coordinates 0i,<Pi , i=l,..,N, and FqIs the 

intensity of polar vortices. 

The stream function is as follows 

y, = -nR2coS0 + M\±^ + J-Zr.\R-. , (3) 

2;r l-cosfc* 2^ i=[ ' 1-cosMj-(^,^0) 
where cosw,. =cos^cos^,. -l-sin^sin^,. cos(^-^,.). The first term in the right-hand side of (3) 

^0 1 -l" c os^ 

corresponds to the solid body fluid rotation. In (3) it is taken into account that y/Q = 2^''^ j — cos0 

1 d ^V^n 

a solution of the equation - . ^ ,^ sin^ , = F n.( 5(0) -5(0-71)) . Contrary to the function 

sm6' 00 o0 ^ ' ^ 
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r 1 

~ 1 — co^O ViSQA in [1] (corresponding to the elementary vortex object in [1]), function Wq is 

an exact stationary solution of the three-dimensional hydrodynamic equations of ideal 
incompressible fluid, when in the absence of the radial fluxes with Vr=Q the radial pressure gradients 

may be not equal to zero ( — ^^0), resulting in that !Pg does not satisfy to these equations. The 

dr 

function Wq corresponds to the polar system for the point CVPs and after the replacement (from 
symmetry considerations) of cos^ by cosm, defines the only elementary 8-singular vortex object on a 
sphere allowed by the hydrodynamic equations. 

It may be shown that (2), (3) correspond to the exact weak solution of the equation (1) when 
di(t) and (pi(t) in (2), (3) agree to the 2iV-dimensional system 

R^sm0i TA^k-A l-cos^M.^ 

(4) 

Here we use !f^from (3), = const for an arbitrary Q and Fq, and cosw^^ = cos for ^ = ^^and 

(p = (Pf^. To prove this (as for the vortex dipoles in [9]), it is necessary to substitute (2) in (1), then the 

result of differentiation to multiply by an arbitrary finite fiinction O(0,(p) and to integrate it over the 
sphere surface. After integration by parts (taking into account the non-divergence of velocity field), 
one gets the system (4) by equating to zero expressions at 0(0, cp) and its same derivatives based on 
the arbitrariness of O. 

Obtained in (4) the exact weak solution of the hydrodynamic equations (1) takes into account 
exactly the sphere rotation effect qualitatively in a different way in comparison to the approximate 
modeling of the background and local vorticity in [1,2,6], particularly. In other aspects equations (4) 
(deduced first time herein fi-om hydrodynamic equations) with the accuracy up to the constant factor 
in the right-hand side of the equation coincide with proposed in [1] from kinematic considerations 
and widely used (see [10]) equations written for any CVP number N>2 composed from point 
vortices. Note that accounting in [11,12] of the sphere rotation effect in the exact solution of the 
absolute vortex conservation equation (1) practically coincides'^ with given in (3). An important 
distinction fi-om [12] is in the use in (3) of the fimdamental solution for the spherical part of Laplace 
operator instead of the operator's eigen functions (as in [12]) which have singularities corresponding 



There is correspondence between first members in right-hand side of (3) and in (4) from [12]. 
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to the point vortices on a sphere. 

It was suggested in [12] a description excluding these singularities while in (4), vise versa, finite 
dimensional dynamics of such 5-singuIar vortices is considered. This difference in the approaches 
corresponds to the manifestation of the "vortex-wave" dualism for an object with continuous (wave) 
and discrete (vortex) properties at the same time. 

Important distinction of (4) fi-om systems considered in [1,2] is the sign minus at Q that just 
corresponds to the exact weak solution (4) of the equation (1) for co from (2). Used in (4) stream 
function form (3) allows exact accounting of the sphere rotation effect in the conservation law (1) and 
corresponds to the general description based on (4) for time evolution of arbitrary systems of point 
(logarithmic) vortices on a rotating sphere. It also allows study of stationary vortex modes without 
restrictions and artificial complications due to conventional (see [13]) stationary reduction of 
equation (1) leading to the necessity of screened (Bessel) point vortices introduction. 



2. The system (4) for N=2 has relative stationary mode with = and — = , e.g. in 

dt dt 



(5) 



case of equal longitudes (p^if) = (P2(t) = (pit) , if the following equality holds 

r'Q(sin^ ^20 ~sin^ O^q) sin^j^Q + sin 
sin^j^Qsin^20 ^^^(^20" ^10-^ 

where 0^ = 6^^ = const , 0^ - 6^^ = const , 0^^ ^ 0^^ . It is possible that each vortex moves along the 
latitudinal circle {0 is the latitude complement) with constant angular velocity 

ItjR^ sin^6'jQ sin^6'2Q 2;rf?^ sin^^Qsin^20^"^(^20 "^10-^ 
Absolute stationary mode with co=Q in (6) may exist only when /(^O or Qt^^O, and relative 

^^ sin^^rv 

with mi^) is possible according to (5) for = - . /■^ and under r'o=0, Q=0, when from (6), one 

gets <y = - — ^ ^ . In particular, for r2<0 and 6^ <9^ one gets fy<0, that 

nR'^ sin^jQsin(^20~^lo) 

corresponds to the widely spread splitting blocking [5] with an anticyclone (of intensity i''^) to the 
north of the cyclone (with intensity Fi) and with intensity of the anticyclone greater than that of the 
cyclone. In that case all the system moves with the angular velocity \(X)\ to the west. It is worth to note 
that in [14] for topographical vortices (in a rotating vessel with the parabolic bottom), it was stated 
similar relationship for intensities of two vortices and the distance from the axis of rotation. In [13] 
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2 

r, sin 0^^ 

also very similar ( = - ^ — in [13]) condition for the relative stationary mode with coi^O was 
^2 sin 

obtained for two vortices on a rotating sphere with the use of Bessel vortices and initially stationary 
reduction of equations (1). 



3. Condition co=0 in (6) for Q^O and /o=0 in (5), (6) defines a new stable absolute stationary mode 
not considered earlier ([13,15], in particular). According to observations [3] it is displayed for 
coupled ACA above the NH oceans, for instance. It is essential that for this stationary solution, less 
intensive (compared to the cyclone) anticyclone with / 2<0 (under Q>0) is always positioned closer to 
the equator than the cyclone. Actually, according to (5), (6) for /"o=0 and co=0 with ^lo S^ts: 

/"j = -;z^?^Qsin 0^ q sin(^20 ~ ^1 0^ ' ^^^^ 

For the typical ACA positions in the NH one can estimate from (6a) {Fj |«10 km /day. It 
corresponds to the existence of an absolute stationary solution without considering the possible polar 
vortices effect, in principle allowing decreasing of that estimate to the usually used empirical value 

10 km I day (see [4]). For the stated new absolute stationary solution, there is not only noted 
qualitative but also quantitative agreement for estimates of the intensities ratio F^iF^ based on 
observations [3] and theoretical estimates according to the equality F^lF^ =-sin^20 /sin^^Q under 
^20 ^^10 ^^"^^ A <0). For example, in comparison with the data (for 1949-2002, in particular) for 

the North Atlantic pair Icelandic cyclone - Azores anticyclone, the mean deviation of the theoretical 
and empirical estimates of the intensities ratio F^iF^, is below 5%. However, respective deviation 
for the Pacific ocean vortex pair Aleutian cyclone - Hawaiian anticyclone is already remarkably 
larger (46%). 

4. Let now Fq , Q^, and sin B^^ ^ sin 6^^ . In that case, from (5) and (6) for a)=0 negative 
F, 

<0(i.e.. corresponds to an anticyclone) is possible only under condition 



F / 

y^= Yj-, <— sin^jo/sin^2o- Then, for example, for ^20^^10' follows from (5) that 

/ 2 

F / 

7q = > (i.e. the polar in the NH is also an anticyclone, see Fig. \)\f -sin ^20 ^ ^10 • 

Relative stationary solution (5) and absolute stationary solution (5), (6) (for cy=0 in (6)) now 
(contrary to the case of Fq = ) may be exponentially unstable with respect to the arbitrarily small 
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magnitude disturbances if inequality (A. 6) (see Appendix) is violated, corresponding to the case when 
/?o = in (A.6). Figure 1 gives examples of the use of the pointed out inequality (for Rq =0 in (A.6)) 
in comparison to the observed data for AC As in the Northern Hemisphere [3]. In particular, if the 
inequality (A.6) (for i?o = ) is met then the linear stability for y « 25 " (Icelandic cyclone) and z « 55 

° (Azores anticyclone) is provided for negative and 1,35<|;;/J<5,24 when according to (5) 

ro«o,3(i,94-|ri|). 

5. In the case of sin 0^^ = sin O^o ^r -T^O, particularly for 02o~^~ ^lo ' ^ more simple inequality 
1 + 4-^ cos <0 is , instead of (7), a condition of the exponential instability of that stationary 

solution with = -1, = • 

To compare this condition with the conclusions on the stability in [6] (based on approximate 
taking into account of the sphere rotation), let's consider an introduced in [6] parameter 

Mr) r, Fr. — Fn. 2 cos 6*1 n 

G = Jr = -j^ c os6>i ^ + -j^ c os6>9^ = const with M const . In particular, -j^ = ^ for 

^0 ^0 ^0 M ^ 

^20= -^"^10 ^iid Yy=-\. For 9^^^<7il2 from the exponential instability condition 

Fr. 

1 + 4-^ cos 6*^0 <0 it follows that the linear instability takes place only when G<0 and |G|<Gq, 

where Gq =8cos^^jq. For instance, Gq is about 0.08 for the pair of vortices near equator with 

« 84.3° according to [6]. Note that in [6] there was obtained quantitatively close oscillatory 
instability threshold value with Gq =0.1. According to (6), the motion of the vortex pair to the west 

with (y<0 (simulating blocking) may take place for Fj > when G < and | G |< 4cos^ 0^^^ , and for 

Fj < under all G > (for G < when | G |> 4cos^ 9^^ ) unlike to [6] (where it is possible only for 

G<0). Thus, conclusions on the stability of relative stationary mode with (o<0 important for the 
blocking realization are found to be essentially dependent on the accounting of the finite intensity of 
polar vortex with Fq ;^ . It would be interesting to compare theoretical conditions for the stable 

blockings possibility obtained here to the results of analysis for polar vortices intensity from 
observations. 
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6. The exact weak solution of hydrodynamic equations on a sphere in the form (4) allows 
arbitrariness in the choice of the polar vortices intensity r,, , including also possibility of evolution of 

Tq in time (T^ (t)) . This arbitrariness may be cancelled if in addition to variability in time for Tq 

variability in time for the sphere rotation frequency is allowed (Q = Q(^)). The requirement of 

preservation of surface integral (for the sphere) of energy E (see (A.l) in Appendix) and angular 
momentum M^(see (A.4) taking into account that in the spherical coordinate system =M^ =0), 

the system of equations (A.5) for variables Q.(t) and (t) can be written out. Accounting of (A.5), 
in particular, leads to the appearance in the inequality (A. 6) of the term Rq j^O, which substantially 
affects on the definition of the stability domains for the stationary mode (5), (6) (when a>=Oin (6)). 
Figure 2a shows stability domain (shown by dark-grey) for the stationary mode (5), (6) when =-l 
. Figure 2b (illustrating sensitivity of the stability domain to the /i variation) shows the stability 
domain when = -1,001 . The same time, without taking into account of (A.5), i.e. for =0 in 
(A.6), the mode with = -1 is unstable for all ^10,^20 when sin^j^ ^ sin ^20 • 

The case with sin 0^^ = sin is considered in section 5 and, in more details, in Appendix. 

7. In the conducted above analysis of the point CVP dynamics on a rotating sphere, non-uniformity 
of the subjacent surface, non-adiabatic processes and other factors (including finite size of the vortex 
core) characteristic for a real atmospheric vortex system were not taken into account. However, it is 
possible to estimate impact of the purely dynamic component of the vortex interactions on the 
evolution of large-scale cyclonic and anticyclonic ACA playing an important role in forming of the 
planetary climate system and its regional pecuharities [3]. Attempts to describe ACA dynamics with 
the help of the point vortices on a rotating sphere taking into account polar vortexes was undertaken 
in particular in [4] but without any assessment of these factors contribution in comparison to 
observations. 

Based on the proposed here theoretical approach the appropriate analysis for the vortex modes 
of ACA pairs from long-term observations [3] was conducted. Analysis of pressure fields allows to 
estimate parameter (ratio of the circulation intensities for a cyclonic and anticyclonic vortices in 

the respective ACA pair) used in the theory. For this purpose, data on the pressure anomaly in the 
each ACA center with respect to its periphery with a given latitude-longitude ACA coordinates (see 
[3]) were used. 

Figure 1 shows examples of the stability domains (indicated by dark-gray color) according to 
(A.6) (for Rq=0) depending on the complement of the latitude for the cyclonic and anticyclonic 
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vortices of different intensity. Cyclonic and anticyclonic vortices intensity ratios specified on Fig. 

1 were realized (see [3]) for the North Atlantic ACAs in winter 1981 (Fig. la) and in winter 1988 
(Fig. lb). Respective positions for Icelandic cyclonic and Azores anticyclonic ACAs are noted on 
Fig. 1 a,b by crosses with the following coordinates: a) 0^ = 17.5°, 47.5°, b) 0^ = 30°, 02 = 
57.5°. According to (A.6) (for Rq=0) and Fig. la the position for the vortex pair of Icelandic and 

Azores ACAs in winter 1981 corresponds to the instable mode, whereas, in winter 1988 on Fig. lb, to 
the stable mode. 

To estimate characteristics of variability and instability of the vortex pairs mode from 
observations the variability of the distance between ACAs along the longitude A9 for the pairs of 
Atlantic (Icelandic cyclone and Azores anticyclone) and Pacific (Aleutian cyclone and Hawaiian 
anticyclone) ACAs was analyzed, in particular. The excess of the standard deviation of Acp relative to 
the mean value over the analyzed period (1949-2002) was taken as a criterion of anomalous 
variability (and instability) for mutual ACAs positions, particularly in winter. Theoretical estimates of 
the stability/instability of a vortex pair mode correspond to empirical estimates of the longitudinal 
(Acp) normality/abnormality with probability 78% for the Icelandic cyclone - Azores anticyclone 
comply with the empirical estimates of the longitudinal (Acp) and with remarkably less probability of 
54% for the Aleutian cyclone - Hawaiian anticyclone. 

One of the factors not included in the proposed above theoretical analysis which must affect 
on the relative dynamics and stability of the ACA vortex system is the temperature gradient between 
the oceans and continents. In relation with this, there was conducted analysis of the observed ACA 
mutual longitudinal positions (Acp) in the NH depending on the temperature difference AT near ocean 
and land in particular for the winter season according to observations (see 
http://www.cru.uea.ac.uk./cru/data/). For the North Atlantic ACA pair the probability of concurrent 
normality and abnormality for Acp and AT was estimated about 67% and about 70% for the North 
Pacific ACA pair for the period 1949-2002. 

According to obtained results contributions of dynamic and thermal factors to the formation of 
anomalous mutual ACA position are comparable, in winter particularly. 

Thus, it was found existence of the stable stationary solution for the Hamiltonian djoiamic 
system of two point vortices of different signs on different latitudes of a rotating sphere in general 
agreement with observations. It is shown that such a mode realization is possible only in the case 
when the more intensive cyclonic vortex is on the higher latitude than the anti-cyclonic vortex. 
Exponential instability criterion for the stationary vortex solution taking into account the polar 
vortices effect was obtained. Correspondence of the theory with the observed and reanalysis data for 
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the coupled quasi-stationary systems of cyclonic and anticyclonic ACAs in the Northern Hemisphere 
over oceans was considered. 

The authors are grateful for useful discussions to G.S. Golitsyn, V.P. Goncharov, V.M. 
Gryanik, M.V. Kurganskiy, G.M. Reznik, V.Yu. Tseitlin, O.G. Chhetiani. This study was supported 
by RFBR and RAS programmes. 



Appendix 

1. Dynamic 2N - dimensional system (4) provides conservation of integral invariants of 
kinetic energy E , angular momentum M and momentum P , per mass unit. Here upper bar denotes 

integration of corresponding values over the sphere surface, and E = ^\^ , M = [rxV], P = V, 

dr 

V = — , r - is a radius-vector in the Cartesian system of coordinates {x,y,z) origin of which is 

placed in the center of the sphere. Rotation of the sphere with angular velocity Q is performed 
around axisz . In the result of averaging one gets P = 0, whereas E and vector M components are as 
follows 

^ ^ r,r,, . i+cosM, 



E=—y y li^in — ^-^-^ + 

8;z- ,.=1 R 1-cosw,.^ 

(A. 1) 



KR'j^, ' 1-cos^, 3 I ° tr ' ' 



+ 

2nR 

= 2Xr, cos^, +4r„ , (A. 2) 

1=1 ~> 

_ N _ N 

= 2^ r,. cos (p^ sin e^ , My^'^Yj ' ^^-^^ 

where 9^ and ^, , i = l,l...N , are the functions of time, defined from solutions of the system (4) with 
corresponding initial conditions. When Q = and = all four values (A. 1) - (A.3) are invariants 
(4), matching with invariants of the dynamic system in [1]. When Q;^0 or ^^0 values (A.3) are 
not invariant, but (A.l) and (A.2) are still invariants of the system (4). 

2. Let's note that in the spherical coordinate system for M and P, we have: M^. =0, 
= RVg, Mg=-RV^; P^=Q , P^^Vg, P^=V^ (here as before r = R,V^=Q). For any Q and F„ 

- M - - 

there is the only independent non-zero integral value Mg (because P^- , = , = ) in 

R 

the form 

_ ,V N 

Mg = TT^R^n - 27^0 - 2j^r.0. - ;r^F. . (A. 4) 

i=l i=l 

Value (A.4) is not invariant of system (4) in the general case, since (A. 4) and (A. 2) are not the same. 

3. Weak solution (4) allows the dependence of Q and F^ on time. From the conditions 

we get in addition to (4) the system 
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2 

fli„ . 2r,-7di'a = —Y^rA, (a.5) 

where a = — X^ln- /-2f^, b = -2M,^-2r, + ——, and M„,=Xr,cos^, is an 

iTrtTiR 1-cos^^ 3 ""T 

invariant of the system (4). The use of system (A. 5) together with (4) eliminates arbitrariness in the 
definition of the value of Fq, allowed when getting a weak solution of equation (1) as a system (4) 

(since in that case we get only conclusion on constancy of the values T. i = 1, 2...N for arbitrary Tq ). 

When considering stability of the stationary modes (5), (6) (for &) = in (6)) taking into account the 
system (A.5), it is found out that stability to extremely small disturbances takes place when the 
following inequality is satisfied 

,2 ^ , ., D , A.. /- -:„4 , „• 4 



Y^A + y^B + C 47o(cosjsin z + ^j sin jcosz) 



r^e 



sin^(z -3;)sin^ zsin^ >> sin^>'sin'^z 



+ Ro, (A.6) 



A = sin^ y(siny-2sinzcos(z-j;)),S = -2sin^ zsin^ >'cos(z-jv), 
C = sin^ z(sinz-2sinjKCOs(z-jK)) 



^0 = 



4/j (sin^ z - sin^ j;)(sinz - sin y) 

;r sin ^ z sin ^ j;(2 + i?j ) 

, 1 + cosy , 1 + cosz ^ 

r,]a ^ + ln^ 2/2 

l-cosj; 1-cosz 

2 

4(-r2 -7o -cosz-/iCOSj) 

(;Ki sin J + sinz)sin jsinz /^(sin^ j + sin^ z) ^ fjSiny-sinz 
2sin(z->')(sin^ z-sin^ jv) ' ^ sin^>'sin^z 2sin(z->')sinzsin>' 

r r 

In (A.6) as previously = —,7o = ~' but = 27iR^Q./T2 . 

Fj 

4. For the stationary mode with c 6*20 = ^ - ^10 h = <P20, having place for ^1 = -1 only, equations 
of the small disturbances are as follows: 



4sin^ipCOS 
^i(l-47ocos6'io) 



sin 

where q)-q)^-q)2 and 6^ are extremely small disturbances of the stationary mode , ^20 = ^ - ^10 > 

^10 ~ ^20 ' the same time corresponding disturbance of the state ^20 the same as for 

— TiR^Q. . 2 1 

according to the invariant M (see (A.2)). In (A. 7), the value = sin 6'jo H is 

F2 4cos^io 

defined from the condition of absoluteness of the stationary mode , O20 - ^~ ^10 5 ^10 - ^20 ■> 
excluding transport of the vortex pair along the latitude, i.e. for = Oin (6). For (A.7), it follows that 
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absolute stationary mode is stable when nR^Qsii ^6^^ IT^ < 0, i.e. when the vortex with intensity 
is anti-cyclonic, and the one with intensity Fj is cyclonic. It means that namely when accounting 
finiteness of the intensity of the polar vortices y^^O , it is possible to have stability of the 
configuration of two cyclonic vortices placed symmetrically relative to equator. The same time, if to 
consider not an absolute but a relative stationary mode (with co^Qm (6)), then for arbitrary with 

f Q < ^- the systemcHCTCMa (A.7) yields exponentially stable solution for pairs of anti-cyclonic 

4cOS^m 

vortices placed symmetrically relative to equator. 

For sin = sin 6*20 , another stationary mode is possible with sin6'jo = sin 6*20 , when 6'jq = 6*20 
and (p^Q - (P2Q = ;z" . In that case, already with necessity, Y^=l . Stability of such a stationary state of 
the pair of identical vortices on the same latitude is defined by a system of equations coinciding with 
(A.7) but with the change of signs in the write hand of (A.7) on inverse and replacing of y^ by -y^, 

where now = ^ ^ sin^ 6^^ ^- . Thus far, stable absolute stationary mode is found to be 

Fj 4cos6'io 

possible for the case of anti-cyclonic vortices only as for y(^=0 , and for Yq ^0 also. The same time, 

for two identical vortices of any sign, it is possible to have relative stationary mode when the vortices 
move with the constant speed along the latitude (i.e. for co^O in (6)). 
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Fig. 1. Theoretical domains of stability (shown by dark-grey) of the vortex pair on a 
sphere according to (A.6) (for R^^ =0) depending on the complement to the latitude for 

anticyclonic (horizontal axis Z) and cyclonic (vertical axis Y) vortices with different 
values of ratio of their intensities : a) =-1.10 and b) =-1.87 . Observed 

positions of Icelandic cyclonic and Azores anticyclonic ACAs are shown by crosses: a) 
for the winter of 1981 with /^ =-1.10; the cross coordinates are y = 17.5° , z = 
47.5°, b) for the winter of 1988 with y^ = -1.87 ; the cross coordinates are y = 30°, z 

= 57.5°. Shaded domains correspond to positive values of /r, = - ratio of the polar 
vortex intensity Fq to the anticyclonic vortex intensity F2 . 
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Fig. 2. Stability domains (shown by dark-grey) for the stationary mode according to (A.6) with 
Y,=-\ (a) and =-1.001 (b). 
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Delta-singular vortex dynamics on a rotating sphere and stability of coupled atmosphere action 

centers 



I.I. Mokhov, S.G. Chefranov, A.G. Chefranov 
Abstract 

Existence of a stationary mode for a Hamiltonian dynamic system of two point vortexes with 
different signs on different latitudes of a uniform rotating sphere complying with observed data is 
stated. It is shown that such mode realization is possible only in the case when the more intensive 
cyclonic vortex has greater latitude than that of the anticyclonic vortex. A criterion of exponential 
instability of the stationary vortex mode taken into account impact of the polar vortexes is obtained. 
Compliance of the theory to observed data and reanalysis for coupled quasi-stationary systems of 
cyclonic and anticyclonic atmosphere action centers above oceans in the Northern Hemisphere is 
considered. 
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